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Abstract 

In the past decade or so observations of supcrnovae, Large Scale Structures (LSS), and the 
Cosmic Microwave Background (CMB) have confirmed the presence of what is called dark energy 
- the cause of accelerating expansion of the Universe. They have also measured its density as well 
as the value of other cosmological parameters according to the concordance ACDM model with few 
percent uncertainties. Next generation of surveys should allow to constrain this model with better 
precisions, or distinguish between a ACDM and alternative models such as modified gravity and 
\ (interacting)-quintessence models. In this work we parametrize both homogeneous and anisotropic 

fSJ ■ components of matter density in the context of interacting dark energy models with the goal of 

discriminating between f{R) modified gravity and its generalizations, and interacting dark energy 
^ ■ models, for which we also propose a phenomenological description of energy-momentum conserva- 

tion equations inspired by particle physics. It is based on the fact that the simplest interactions 
between particles/fields are elastic scattering and decay. The parametrization of growth rate pro- 
posed here is nonetheless general and can be used to constrain other interactions. As an example 
■ of applications, we present an order of magnitude estimation of the accuracy of the measurement of 

I these parameters using Euclid and Planck surveys data, and leave a better estimation to a dedicated 

work. 



2 ■ 1 Introduction 



Nowadays it is a well established fact that according to the Einstein theory of gravity ~ 73% of the 
mass and energy in the Universe is in a strange form with unusual properties inconsistent with any 
^ ■ type of matter known to us. It is generally called dark energy (see e.g. jTl[2l[3ll3] for recent reviews). 
. In the last two decades or so numerous models have been suggested to explain this mysterious and 
[ dominant constituent of the Universe. The majority of these theories can be classified in one the 

■ following three categories: 1) Models based on a scalar field e.g. quintessence [51 El [3 |8] and its 

■ variants such as K-essence O [TUl [TT| [T2| [T3] in which the kinetic term in the Lagrangian has a non- 
standard form, and varying neutrino mass models |14^ [T5] in which the accelerating expansion of the 
Universe is generated by the variation of neutrinos mass due to their interaction with a light scalar 
field; 2) Modified gravity models in which dark energy is explained as the deviation of gravitational 
interaction from Einstein theory of gravity. Examples of such models include scalar-tensor [16\ [T71 [T8] 

X : and f{R) gravity [ISl EQl EZ] , Chameleon [SH [22l [23] , and DGP [S]; 3) A cosmological constant - 

■ introduced by Einstein himself [25] and interpreted by Lemaitre as the energy density of vacuum [26] . 
It is phenomenologically the simplest of three categories, and is still the best fit to all observational 
data |65 [ 145 1 166]. However, naive estimations of vacuum energy are ~ 42 to 123 orders of magnitude 
larger than the observed dark energy |27j . For this reason, alternative explanations have been explored 
even before the observation of the accelerating expansion of the Universe [5l [6]. The main task of 
cosmologists today is discriminating between these models, in particular distinguishing the first two 
categories mentioned above from a cosmological constant. 

A notable difference between a cosmological constant and some of alternative models is the presence 
of a weak interaction between matter and dark energy. Pure quintessence models, in which there is no 
interaction between the scalar field and matter are somehow pathological because all known fundamen- 
tal particles, including neutrinos which have very small couplings, interact non-gravitationally with 
some other particles. Even axiomatic weakly coupled particles such as axions [281 [29] are expected to 
interact with gauge bosons such as gluons. Fields in candidate extensions of the Standard Model are 
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related to each others by symmetries, thus either by gauge interaction or by mass mixing. On the other 
hand, if a field such as quintessence interacts only with gravity, then naturally it should be considered 
to belong to the gravity sector. An example of such fields is dilaton which was first introduced in the 
context of Kaluza-Klein model for the unification of gravity and electromagnetic forces \30\ [3T] , and 
is also associated with conformal gravity models, see e.g. [32j and references therein. But gravity is 
a universal force and interacts with all other particles. Thus, in contradiction with the assumption 
above, the quintessence field must have an interaction with other particles. In fact dilaton does have 
non-minimal interaction with other species, see for instance [331 134j . This makes the task of finding a 
candidate for a non- interacting quintessence field very difficult. A more problematic issue with pure 
quintessence models is the fact that they do not solve the coincidence problem of dark energy, i.e. why 
it becomes dominant only at late times and after galaxy formation. Interacting quintessence models 
in which the quintessence field has a weak interaction with some matter species, in particular with 
dark matter, can solve or at least soften the huge fine-tuning of dark energy density with respect to 
matter in the early Universe [35\ [36] . 

In modified gravity models the deviation from the Einstein theory of gravity can be, either written 
explicitly, or presented by introducing new fields, usually scalars, in the matter sector. The first 
presentation is called Jourdan frame and the second Einstein frame. Because in the latter case the 
model looks very similar to an interacting quintessence mode0, it is necessary to find a proper definition 
that discriminates between what is called modified gravity and what is called interacting quintessence. 
In modified gravity, the scalar field is usually related to a dilaton field, thus it has a geometrical origin 
and arises from a broken conformal symmetry [16j . For this reason, the scalar field always interacts 
with the trace of energy- momentum tensor of matter |67j . The situation is not so straightforward 
for interaction between matter and scalar field in interacting quintessence models for which various 
types of interactions are considered in the literature, see for instance [70j . In many of these models 
in analogy with modified gravity, in particular i{R) models, the interaction term is considered to be 
proportional to the trace of the energy-momentum tensor of matteiH. 

In this work we try to determine the interaction between dark matter and dark energy in a collisional 
description of interactions inspired by particle physics. Using the Boltzmann equation with a collisional 
term and some results from studies of the microphysics of dark energy condensate [7T], we show that 
the interaction can be described only approximately by spacetime dependent functions, and in general 
one needs the distribution in the phase space f{x,p) where p is the 4-momentum, see Sec. [3] for 
more details. However, at present and for foreseeable future, we cannot observe the phase space 
distribution of dark energy. Considering this fact, we use thermodynamical description of average 
energy-momentum and velocity to obtain approximate covariant expressions for interactions between 
matter and a scalar field as dark energy. This leads to a modification of energy-momentum conservation 
equation which explicitly deviates from modified gravity. Their difference can be used as a mean for 
classifying models and discriminating these two categories in the data. 

On the observational side, one has to find the best way of parametrizing cosmological evolution equa- 
tions such that they admit discrimination between at least the three major categories of models 
discussed above. In preference they should not depend on the details which are neither well under- 
stood nor can be targeted with the precision of present and near future surveys. Observations show 
that dark energy has negligible clustering (see e.g. [72l [73l [TH [75] for latest results). Therefore, its 
dominant contribution is in the homogeneous component of the Einstein and conservation equations. 
It also affects the evolution of anisotropies mainly through their dependence on the background cos- 
mology. For this reason, irrespective of the way we measure the equation of state of dark energy - 
from observations of supernovae that are only sensitive to background cosmology, or from observations 

'^Note that when we talk about interacting quintessence models we mean models in which the scalar field interact 
with some other constituents of the Universe. All quintessence models have a self-interaction which is not explicitly 
considered in the formulation presented in this work 

^For the reasons described in detail in Sec. [3l when we talk about the interaction term, we mean the modification of 
energy-momentum conservation equation due to an interaction. 
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of matter perturbations by measuring lensing or galaxy distribution - we must extract parameters of 
background cosmology to determine the contribution of dark energy. Consequently, it is crucial to un- 
derstand how different models affect this component through a proper parametrization that facilitates 
the discrimination between various models. This is another goal of the present work. Although there 
are few popular parametrizations \76\ \77\ \7E[ [79] in the literature, specially for testing modification 
of the Einstein theory of gravity at large scales, as we will show in this work, they are not suitable 
for discriminating between modified gravity and (interacting)-quintessence models. We should remind 
that for ACDM model the growth rate f is roughly scale independent. Therefore, observation of 
the violation of this property would be a clear signature of inconsistency with standard cosmology. 
But the measurement of f and the expansion rate H by themselves are not enough for discriminat- 
ing between modified gravity, quintessence, and interacting quintessence, and a parametrization that 
does not depend on the details of these models is necessary to highlight their differences. Evidently, 
one can simply fit the data with models and compare their goodness of fit. But, this does not take 
into account the degeneracies and similarities. Therefore, a smart parametrization and better data 
analysing methods are necessary. Moreover, the fact that most popular modified gravity models can 
be formulated as an scalar field theory means that their differences from (interacting)-quintessence 
must be understood and the parametrization must be performed in a way that it highlights these 
differences and help discrimination. 

In this work we propose a new set of parameters to describe, in a model independent way, the effect of 
an interacting dark energy on the evolution of the expansion rate of the Universe and another set of 
parameters for the growth rate. These quantities are the most sensitive measurables for discriminating 
between dark energy models. Consequently, the ultimate goal of various measurement methods is 
to constrain cosmological and dark energy models by measuring one or both these quantities. For 
instance, galaxy distribution and lensing surveys determine the power spectrum of fiuctuations for 
one or multiple redshift bins. Future large surface and sensitive spectroscopic surveys such as Euclid 
allow to determine the matter power spectrum for a statistically significant number of redshift bands, 
and thereby extract the growth rate, see e.g. |80j for the methodology applied to The WiggleZ Dark 
Energy Survey. The BAO measurements determine the expansion rate and angular diameter distance 
at one or multiple redshift bins, and can thereby estimate variation of these quantities. Supernovae 
data measure the expansion rate directly, and the variation of H{z) can be extracted. Therefore, 
parametrizations that we will discuss here are relevant for all measurement methods in cosmology. 

In Sec. [2] we present a newparametrization for Friedman equation in the context of a general in- 
teracting dark energy modelj. In Ref. [37] we defined a quantity B{z) oc dp/dz and proposed it for 
the measurement of the equation of state of dark energy defined as u; = Pde/Pde where Pde and pde 
are pressure and energy density of dark energy respectively. It is specially suitable for measuring the 
deviation from a cosmological constant, see Appendix [A] for the definition of B{z) and a review of 
its properties. In addition, we argued that in what concerns the sign of 7(2) (see equation ([2|) below 
for its definition), this quantity has distinct geometrical properties which make it less sensitive to 
uncertainties of other quantities such as Hq or respectively the present value of Hubble constant 
and the density fraction of matter. The sign of 7(2) is the discriminator between what is called phan- 
tom models which have w < —1, and normal scalar fields (quintessence) models and a cosmological 
constant for which w ^ —1. 

Using this parametrization and properties of B(z), we show that in presence of an interaction between 
dark energy and other components, one obtains different estimation for 7^g'^ (see next section for its 
definition) from H{z) and from B{z) when data is analyzed with the null hypothesis of a ACDM 
model as dark energy. In this way one can predict the sensitivity of surveys to interacting dark energy 
models in a model-independent manner. Then, we discuss the properties of parameters for each 
category of models, their differences, and how this information can be used to discriminate between 
various dark energy models. In Sec. [3] we describe phenomenological interactions for interacting 

^After the submission of this paper a similar parametrization for the Friedman equation is reported independently 
by [38]. 
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quintessence models and compare it with modified gravity. This leads to an approximate description 
for non-gravitational interactions between dark matter and dark energy. 

In Sec. [Hwe present evolution equations of over-density and velocity fields in each category of models 
for the interactions obtained in Sec. [3l Then we describe how one can discriminate between interacting 
quintessence and modified gravity models by using matter power spectrum and its evolution, i.e. the 
growth rate of anisotropies. Because the growth rate plays a special role in the discrimination between 
various dark energy models, in Sec. [5] we parametrize its evolution, and as an example of application, 
we obtain an order of magnitude estimate for the discrimination ability of the Euclid mission [l^ 
that measures both parameters of the homogeneous component (the background cosmology) and 
the evolution of growth rate of matter anisotropies. In addition, we compare our parametrization 
with other parametrizations that can be found in the literature which are usually designed to test 
the Einstein theory of gravity. Conclusions and outlines are summarized in Sec. [6l Properties of 
the functions B{z) (and A[z)) are reviewed in Appendix [XI Fisher matrix for dark energy without 
parametrization of its equation of state w{z) is described in Appendix [Bj A summary of covariant 
formulation of a classical scalar fields as a perfect fluid is given in Appendix O In Appendix [D] we 
calculate an approximate analytical solution for the growth rate of matter anisotropies. 

Here we must emphasize that the predictions for future missions obtained in this work are only 
representative and order of magnitude estimations of what is expect from future surveys. They should 
be considered as a QD (Quick and Dirty), hand-shaking predictions. Their purpose is only to show that 
it is possible to measure the new parameters with reasonable uncertainties. A proper prediction for 
future observation projects needs detailed consideration of instrumental response, simulation of data 
analysing procedure, and understanding of the sources of systematic and statistical errors. Fulfilling 
these requests necessitates a dedicated investigation which is out of the scope of the present work that 
targets theoretical issues related to the discrimination between various dark energy models. In fact, 
a number of authors have performed predictions for uncertainties of various measured quantities by 
future missions, see for instance (Ml El [H21 ES] . They usually consider models that can be classified 
as modified gravity according to the classification criteria discussed in Sec. O Nonetheless, some of 
their parameters can be related to quantities defined in this work, thus their predictions can be used 
to obtain a rough estimation for the expected uncertainties for the new parameters. 

Throughout this work we use Einstein frame for modified gravity models unless it is explicitly specified. 
In this way, a unified description can be made for all interacting dark energy models based on a scalar 
field formulation. 

2 Friedman equation in interacting dark matter-dark energy models 

Apriori the measurement of the equation of state of dark energy is simple. It is enough to measure 
the expansion rate of the Universe H{z) = a{z)/a, or a quantity related to it such as the luminosity 
distance D{z) at different redshifts. Then, by modeling known constituents of the Universe as non- 
interacting perfect fluids, one can fit the data and measure the effective equation of state of dark 
energy Weffiz), defined as Pef f{z) / pef f{z). The suffix "eff" is used to remind that pressures and 
densities obtained in this way can be effective quantities rather than physical pressure and density of 
constituents, because we have neglected any interaction between components. Therefore, from now 
on effective quantities mean quantities determined from data by considering a null hypothesis. 

In practice however the life is not so simple. The density of a perfect fluid changes with redshift as 
(1 -|- 2)^'^ (7 is defined in Therefore, at low redshifts when z — )• 0, the total density is not very 
sensitive to the value of 7 or equivalently w{z) and their variation with 2, see Appendix |A] for more 
details. This statement is independent of the type of data or proxy used for determining H[z) or 
D{z). On the other hand, at high redshifts where H{z) is more sensitive to the equation of state, 
dark energy is subdominant. Moreover, it is more difficult to measure H{z) and D{z) at higher 
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redshifts and measurement uncertainties can make the estimation of w{z) and its evolution unusable 
for discrimination between models. 



If constituents of the Universe do not interact with each other, Friedman equation which determines 
the evolution of expansion function a{t) can be written as: 

^ = Pi5l = n^^i + ,f + + z)4 + QKil + zf + + zf^^^\ p,{z) ^ ^ (1) 

Hq pq SttG 

m = cold dark matter, b = baryons, h = hot matter, k = curvature, and de = dark energy 

In this class of models matter and radiation densities evolve only due to the expansion. This is a good 
approximation for all redshifts z < Zcmb ~ 1100. 

In interacting dark energy models matter and radiation terms in the right hand side of the Friedman 
equation ([T]) can contain an additional redshift-dependent factor: 

= = QiTiizMl + zf"/* i = m, b, h,k, and de (3) 

Ho Pco Y 

Without lack of generality we assume that J^de = 1 and all redshift dependent terms are included 
in 7(2). In quintessence models the coefficient of the curvature term also is constant because it is 
assumed to be related to geometry/gravity and independent of the behaviour of other components. 
At present observations are consistent with only gravitational interaction between various components 
in ([3|), thus additional interactions must be very weak. By definition and without lack of generality 
we consider J^i{z = 0) = 1. Observations also show that ilk ~ 0, therefore throughout this work we 
assume fi^ = unless it is explicitly mentioned. Note that in the case of modified gravity models, a 
parametrization similar to ([3]) can be defined both in Einstein and Jordan frames. 

A simple example for which an approximate expression for J-i{z) coefficients can be found is a model 
with a cosmological constant as dark energy and a slowly decaying dark matter. The decay remnants 
are assumed to be visible relativistic particles [39]. In this case: 

|J ^ nrnii + zf exp(^^) + n^ii + zf + + zf + + zf (^i - exp(^^)^ + 

(4) 

J™(i)«exp(^^^) + (l + z)(^l-exp(^5LZl)^, r > tq, = J'h = I l{z) = (5) 

where r is the lifetime of dark matter and tq is the age of the Universe. It is demonstrated that in this 
example, if the decay /interaction of dark matter is not considered in the data analysis, a w^jf < — 1 
can be obtained for dark energy, see [40t H2] for more details about the set up and the proof. 



Note that in ([5]), we have included the contribution of relativistic remnants in J^^n- However, as 
this component has a redshift dependence similar to hot matter, it also makes sense to consider it 
as hot matter and add it to hot component. It is even possible to add this term to dark energy 
contribution, as long as it is small and induces only a slight deviation from a cosmological constant. 
In this case, one can show that the effective dark energy will have w^jf < — 1 [30]. The reason for 
such freedom is that we do not measure or take into account the decay remnants. This example 
clearly shows that parametrization ([3]) is not unique when all the components and their interactions 
are not know. Therefore, one has to be very careful about degeneracies when data are analyzed and 
interpreted. In particular, prior assumptions such as stability of matter and radiation components can 
affect measurements and conclusions. This example also show that for ruling out ACDM model, it is 
enough to prove that at least one of J^j 7^ 1, or jde 7^ cEI- 



^This statement is true if baryon pressure is negligible. Future surveys can be sensitive to small baryon pressure. In 
this case it must be taken into account before any conclusion about ACDM model is made. 
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Extension of this example to quintessence models without coupling to matter is straightforward and one 
simply needs to consider j{z) 7^ 0. A more interesting extension is to assume that the quintessence 
scalar field is one of the remnants of the decay of dark matter, which during cosmological time 
condensates and makes a classical quintessence field. In this case, it has been shown |36l [85} [7T] 
that coefficients J^m, J^h, and equation of state of dark energy w{z) (or equivalently 7(2;)) are not 
independent. However, their relations are too sophisticated and cannot be described in an analytical 
form and numerical techniques should be employed |36) . 

According to ([S]): 

Tm{z) >J'm{z = 0) (6) 

and because r » tq, J^i{z) coefficients are close to 1 at all redshifts. In general, for an interaction 
which transfers energy from dark matter to other components, the inequality ([6]) is applied because 
at high redshifts one expects a larger contribution of dark matter in the total density than in a non- 
interacting model. Inversely, if energy is transferred from other components, for instance from dark 
energy, to dark matter: 

Tm{z) <J'm{z = 0) (7) 

An example of such models is scaling dark energy \55\ [56] in which at early times dark energy has 
a much larger contribution in the total energy density, but it gradually decays to dark matter and 
only recently its equation of state approaches w ~ — 1. Another example is the class of models called 
early dark energy. Although the original model |57 [ [55] is a pure quintessence/k-essence, there are 
variants of this model in which, there is an interaction in the dark sector [60] or between dark energy 
and visible sector |61] . 

In models with elastic interaction between two sectors, no energy is transferred between them, and 
J^m{z) = J^h{z) = 1. Nonetheless, the phase space of matter and dark energy in these models can 
change and thereby w^e can depend on z. 

For f(i?) modified gravity models homogeneous Einstein equations and energy conservation equation 
in Jordan frame are [67] j^l 

i 

i 

Ir 

Pi + 3Hpi = - {pi - 3Pi), i = m, h, k 

where Dot and prime mean derivation with respect to comoving and conformal time, 

respectively. Subscript R means derivation with respect to scalar curvature R = R^ud^^ ■ We remind 
that at linear order the effect of matter perturbations on R is zero, thus R only depends on z and 
the effect of f(i?) 7^ on the evolution of perturbations manifests itself by changing the background 
cosmology. 

""When equations apply to both dark matter and baryons, we indicate them collectively with subscript m. 
^Here we have written Einstein and conservation equations in Jordan frame because they lead to expressions for Ti 
coefficients which are explicitly very different from quintessence case. 



f" 

la. 

2a2 



(8) 
(9) 
(10) 
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After solving density conservation equation (fTO|) , Friedman equation ([8]) can be written as the following: 



(11) 
(12) 

(13) 

I, -. = 4 (14) 

Equation ([13]) is the energy density of effective dark energy in f(i?) gravity models. Similar to 
quintessence models we can assume J-de = 1- The only explicit difference between ()12p and the 
same equation for an interacting quintessence model is the presence of a nontrivial coefficient for the 
curvature term if ^ 0. Nonetheless, the evolution of coefficients J-i{z) with redshift is different 
from their counterparts in interacting quintessence models, in particular from models in which en- 
ergy is transferred to dark energy at low redshifts, see equation 1^. In fact, the function f(ii) is not 
completely arbitrary and must satisfy a number of constraints. Notably, f(i?)||^|^o — to make the 
model consistent with Einstein theory of gravity in mild or strong gravity fields, and f /j > due to 
stability constraint [Mj. Under these conditions: 

Ti{z) > Ti{z = 0) (15) 

Comparing (fT5]) . Q, and ([7]) one can immediately conclude that the measurement of Tm{z) and its 
evolution with redshift can discriminate between dark energy models in which energy is transferred 
from dark energy to dark matter such as scaling models, and in modified gravity models. But it 
cannot discriminate modified gravity from models in which energy is transferred from dark matter to 
dark energy such as the model discussed in [Ml EHl El] . To discriminate the latter and other models 
of this category from f/j modified gravity, the coefficient of relativistic (hot) component J-h{z) and its 
evolution must be measured. Evidently, such measurements are very difficult. For instance, one has to 
measure very precisely the temperature of CMB at high redshifts or H{z) at a large number of redshift 
bins and fit the data with 7^ 1. In Einstein frame the evolution of matter density is the same as 
in equation (jlip [67], but evolution equation of hot matter is similar to ACDM. In what concerns the 
discrimination from interacting quintessence what is discussed from Jourdan is applicable. 

2.1 Model-independent discrimination of interacting dark energy models 

In this section we show that if ACDM or a simple quintessence are considered as null hypothesis, 
measurements of effective dark energy density and effective equation of state from H{z) and the 
function A{z) defined in Appendix [A] separately, give different values for these quantities if dark 
energy interacts with matter. Similarity of J-'miz), specially if the curvature of the Universe is zero, 
means that we cannot distinguish between interacting quintessence and modified gravity models in a 
model-independent manner - except for the cases explained above. For this reason in this section we 
only study the discrimination between interacting dark energy models parametrized as in equation 
(|12p and a cosmological constant and/or non- interacting quintessence. 

For analyzing cosmological data, ACDM with a stable and non-interacting dark matter is usually 
used as null hypothesis. Nonetheless, the methodology explained below is not sensitive to redshift 
dependence of 7^6; and we can consider the more general case of non-interacting quintessence as 
the null hypothesis. The expansion of the Universe for such cosmologies is ruled by equation ([T]). 
Therefore, we rearrange terms in equation ()12p such that it looks similar to equation ([T|). Then, we 



J'i(z) 





= Pi{z = 




Pc{z) 


Hi 


PcO 




3 


Pde 


SttG ] 



1 + f/j(z = 0) 



1 + Uz) 



iR 6 a-^ 



l + ffl(i?(z = 0)) 
l + fR{R{z)) 



0, Wh = ■ 
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determine effective quantities which are measured by fitting a ACDM or a non-interacting quintessence 
model to data: 

^ = Y1 ^^(1 + + E ^^i^ii^) -!)(! + ^f^' + ^de{l + zf^'^'^''^ (16) 

In null hypothesis model only is redshift dependent and 7^, i = m, h, k are constant. By comparing 
with ([1]) the effective contribution of dark energy is expressed as: 

f](f](l + zf^'"f^'^ = J2 ^ii^ii^) - + ^f^' + ^^.^e(l + zf-^^^^^^ (17) 

i 

In both interacting quintessence and modified gravity models coefficients Fj's are defined such that 
Fi{z = = 1^ therefore at z = the first term in (fT7|l is null and we can separate ^eff and ^de{z)'- 

niff=nde, ^if){z = 0)=-fde{z = 0) (18) 

^S^^) - ^ ^wiTI^ (19) 



where superscript (H) means measured from Hubble constant H. 

Suppose we can also measure A{z) defined in ()109p . and use it to determine the effective density and 
equation of state of dark energy. For an interacting dark energy model parametrized according to (jl6p 
quantities B{z) and A{z) are: 



Biz) = 



1 dp 



3(1 + z)2po dz 

J2 ^^ Uh^) + (1 + Z)^) (1 + ^)'^^^^'^ + ^de{w{z) + 1)(1 + z)3(7.e(.)-l) (20) 



i=m,h,k 



A{z) = Biz)- n^l + z)'^^^-'^ = 

i=m,h,k 

J2 ^^ UiJ'iiz) - 1) + (1 + ^)^) (1 + ^)'^^'~'^ + ^deiwiz) + 1)(1 + z)3(^<^^W-^h) 
i=m,h,k 

Using ()109p in Appendix lAl as the definition of Aiz), we find the following expression for its parameters: 

nif^iwij\iz) + l)il + z)'^'fh^^ = Y^,UiMz)-l) + il + z)^yi + z)'^' + 

i 

n^eiwiz) + 1)(1 + z)37deW = (1 + ^^Aiz) (22) 



where superscript iA) means measured from Aiz). Equations (|17p and ()22p are fundamentally differ 
ent. In particular: 

u;^;;}(z = 0) + 1 



Ol-) _ E.^.^^^+^.e(^^ = 0) + l) 



which in contrast to ^),fj, in general is not equal to ^.^f,. Equality arises only when J-i do not vary 
with redshift i.e. J^j = 1 at all redshifts. This condition is satisfied by the null hypothesis ACDM and 
by non-interacting quintessence models. Therefore, assuming that and are known (e.g. from 
CMB), simultaneous measurements of Hiz) and Aiz) at even one 2; > is apriori enough for testing 
the presence of an interaction between dark matter and dark energy independent of the underlying 
model. Evidently, in practice the measurements must be performed at many redshift bins to improve 
statistics and to compensate for measurement errors. 
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Apriori one can use other quantities such as angular diameter distance Da or luminosity distance 
Di which are easier to measure rather than A{z). However, both these quantities are functional of 
H{z) - through integration of Thus, in general they do not have an analytical expression. 

Besides, their derivatives depend on J^j's only, in contrast to (p2|) that depends on both Fis and 
their derivatives. Therefore, ^eff and 7e// obtained from dDA/dz or dDi/dz will be equal to ones 
determined from H{z) irrespective of the underlying cosmology. This shows that the function A{z) 
(or equivalently B{z)) introduced in [37] has special properties and is well suited for discriminating 
between dark energy models. It can be measured from supernovae data, see |37j for the methodology. 
As for LSS data, one needs to determine both H[z) and its evolution dH{z)/dz to be able to calculate 
A{z), for instance from the BAO and the power spectrum of matter fluctuations |51j . This is not an 
easy task. As an example consider supernovae observations that measure the luminosity distance Dl 
to a supernova from its standardized apparent magnitude. The angular luminosity distance Da is 
related to the luminosity distance, see (jll4p . To determine dDA/dz apriori one can use the measured 
Da, and determine its derivative (slope). However, due to scattering and discreteness of data, such a 
measurement will have large uncertainties. The same problem arises for dH[z) jdz or Ai^z') because they 
depend on derivatives of D^, see equations (jllOp and ()113p . Nonetheless, there are various methods 
such as binning of data, using a fit in place of discrete data, etc. that allow to improve the estimation. 
Near future large area surveys such as Euclid 03], BigBOSS LSST [S3] will be able to determine 
these quantities with relatively good precision, see also Sec. [5] for measurement methodology. In 
particular, large surface spectroscopic and lensing surveys such as Euclid are able to determine the 
variation of total density with redshift dp jdz oc B{z) with good precision. In Appendix [B] we obtain 
the Fisher matrix for dark energy parameters without considering a specific parametrization for the 
equation of state w(£). 



2.2 Discrimination precision 

Measurements of cosmological parameters show that w°^^ ~ — 1 irrespective of which proxy or mea- 
surement method - supernovae, CMB, or LSS has been used. This means that \Ti(^z') — 1| ~ and 
dJ-i{z)/dz w 0. Moreover, addition of J-i{z) to the model increases the number of parameters. Giving 
the fact that we have essentially two observables: H[z) and one of Da{z), Dl{z) or B{z), greater num- 
ber of parameters means also greater degeneracy, thus more uncertainty for discrimination between 
ACDM, a non-interacting quintessence, and interacting dark energy models. 

One way of measuring the presence of interaction without having to fit data to the large number of 
parameters in equations (fT6|) and ([2T]) . is to measure how different ^[ffj, ^l^/' '^e//(-^) ^^^^ 

because as we discussed in the previous section, when J^i ^ 1 these effective quantities are not the 
same. To this end, a natural criteria is: 



^^if](u;^g(z) + l)(l + z)^V/W 



This quantity can be explained explicitly as a function of j , J^j , 7^ , and is zero when Fi = 1, dFi /dz = 
0. Note that we have chosen expression ()22p for comparison rather than ()17p because it is not possible 
to determine ^^g^j in a model independent manner, see equation ([23]) . By contrast f^eJ7 ~ ^de, 

thus %fj and thereby w^j-j can be determined without any reference to F coefficients. In [37] we 
suggested to use the sign and evolution of A{z) to discriminate between dark energy with 7(2;) ^ 
and a cosmological constant. Here Q{z) plays a similar role for discriminating between interacting or 
non-interacting dark energy. 

Assuming that and are determined independently and with very good precision, for instance 
from CMB anisotropies with marginalization over 7^ei © can be determined from the measurement of 
H{z) and B{z). The latter can be measured from whole sky or wide area spectroscopic surveys data 
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such as Euclid, or multi-band photometric surveys such as DES. Evidently determination of B{z) that 
depends on dH/dz is very difficult. However, it is easy to see that there is no other quantity which 
can be measured more easily and discriminates between ACDM and dynamical dark energy models 
with a better precision. For instance, the BAO method determines H{z) and Da{z) directly. But, 
Da{z) depends on w^z) or equivalently 7(z) through an integral, see equation ()114p . Therefore, it is 
less sensitive to the variation of 7(2) with redshift. This is analogue to binning a data. Evidently, a 
binned data is less noisy and has a smaller uncertainty. But, if the goal is to measure the variation of 
data, the binning can completely smear out small variations. Therefore, irrespective of methods and 
measured proxies, we are limited by inherent properties of the physical system. In this respect, the 
precision with which Q{z) can be measured gives the ultimate sensitivity of an observation/data set 
to deviation from ACDM. 



3 Interactions 

In the previous section we used Friedman equation for parametrizing interaction between matter 
and dark energy. Evolution of their densities is ruled by energy-momentum conservation. But, in 
presence of non-gravitational interactions between constituents the energy-momentum tensor of each 
component Tj^'^ is not separately conserved, and conservation equation can be only written for the 
total energy-momentum tensor T^'^ defined as: 

T"'^T.TZree)+TZ (25) 

i 

Ttu=T.T^Zree)-,u + TZ;u = ^ (26) 
i 

where ^^j^j^gg) is the energy-momentum tensor of component i in absence of interaction with other 
components, i.e. T^^^^^^^y^ = 0, and T-^^^ is the energy-momentum tensor of interactiorJll, and T^^^.^ = 0. 
In the literature on interacting dark energy models (see e.g. [35]) when only two constituents - matter 
and dark energy - are considered, the energy-momentum conservation equations are usually written 
as: 

T^;. = Q^ T^% = -Q'^ (27) 

for an interaction current Q^. Comparing ()26p and with ()27p . it is clear that tensors in the left hand 
side of equations in ()27p do not correspond to free energy-momentum tensors, and along with Q^^ 
they are obtained somehow arbitrarily by division of (|26p . In fact, equations in (|27p are inspired by 
perturbation theory in which for each perturbative order, the right hand sides of these equations are 
estimated by using quantities from one perturbative order lower. Thus, they constitute an iterative 
set of equations from zero order (free) model in which = 0, up to higher orders. This approach is 
not suitable for dark energy where we ignore, not only interactions but also the free model. Therefore, 
a more general expression should be used: 

rp^lU _ _f)iJ. rp^U _ _f)IJ. rp^W _ ^fl ^fj. f^o\ 

In these equations matter and dark energy tensors Tm and Tl^^ have the same expression as in the 
absence of interaction, but with respect to fields which are not free. These expressions can be justified 
by considering the Lagrangian of the model. In Einstein frame the Lagrangian for a weakly interacting 
system can be divided to free and interaction parts: 

£ = ^A + A„j (29) 

% 

''Non-gravitational interactions between cosmological constituents must be weak. Therefore, separation of energy- 
momentum tensor to free and interaction component is allowed. 
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Considering only local interactions, in the dynamics equations for the fields partial derivative of Cint 
with respect to each field determines the interaction term. Dynamic equations can be related to 
energy- momentum conservation equations ([27|) [36] . Therefore, interaction currents Qm and are 
generated by partial derivatives of Cint with respect to the corresponding field. 

In the previous section we explained that the scalar field in scalar-tensor modified gravity models is 
related to a dilaton. Consequently, the interaction term is proportional to the trace of matter, see 
equation ()10p for an explicit example of f(i?) models. In this case there is no interaction between 
scalar field and relativistic particles, and it can be shown that Qm = — Q(o [SZ|; i-e. T-'^'^ v ~ ^ ^^"^ 
conservation equations in ()27p can be used. Interaction current for these models can be written 
as: 

Q'^ = C{^)Tmd^''p (30) 

where Tm = g^yTm . In f{R) models the coupling C is a constant. Here we consider (/3-dependence to 
cover more general cases. Some authors have also considered oc TmUm for interacting quintessence 
models |70j . In fact, the interaction current of interacting dark energy models in literature is usually 
considered to be oc pm = 7m for cold dark matter i.e. similar to what is obtained for f{R) modified 
gravity models [67] • However, giving the fact that these models share some important properties with 
modified gravity models, such as the absence of interaction between relativistic matter and scalar field, 
we classify them in the modified gravity category. In fact, interactions in interacting quintessence 
models can be more diverse than this simple case. In the rest of this section we describe how they can 
be formulated without considering their details. 

In the context of quantum field theory, the Lagrangian C can be easily written for various types of 
fields and their interactions, see e.g. [7T]. But these formulations are usually complicated, and are 
necessary if the microphysics of dark energy models is studied. There are various ways to write C 
and/or T^'^ with respect to macroscopic quantities which are apriori measurable from cosmological 
observations. For instance, one can use a fluid description for components. The Lagrangian of a fluid 
is deflned as [83] : 

Cf = ^{P + p)g^,u^u'^ + \{P- p)g,,.g^- + \g^,W'^ (31) 
p = K + y, P = K-V (32) 

where K and V are respectively kinetic and potential energy, and H^^ is the traceless shear tensor. 
Note that if we impose the traceless condition on the Lagrangian, the last term in the right hand side 
of equation (|3ip becomes zero. Therefore, this term must be considered as a Lagrange multiplier, and 
traceless condition is imposed after determination of T^^^ [SI]. It is easy to check that the Lagrangian 
leads to the familiar expression for the energy- momentum tensor of a fluid: 



rjijJLV _ 



Tf = {p + P)u^u'' - g'^'P + (33) 



Transformation of a Lagrangian written with respect to fields to a fiuid description is easy, and one 
can determine the energy-momentum of interaction Tj^^ and the current defined in ()28p directly 
and without any ambiguity, see Appendix O However, their descriptions as a function of density 
and pressure of the fiuid depend on the self-interaction potential V{ip). For instance, a Higgs-like 
interaction between a scalar and a fermion oc ip'ip'ilj is described as oc {p-4) — P^){Pip — PipY^'^ if Vi'f) oc Lp^ , 
and as oc {p^—P^){pip—Pip)^^^ if V{(p) oc (p^. Therefore, when the objective is a general parametrization 
of interactions without considering details of the underlying model, this type of description is not very 
suitable. 

A more serious problem of fluid description of interaction Lagrangian is the fact that conservation 
equations in (|28p are equivalent to field equations and can be obtained from them [36]. Therefore, they 
do not contain quantum processes such as decay and scattering. It is well known that the Boltzmann 
equation plays the role of intermediate between quantum and classical description of interacting sys- 
tems [HU [871 [881 Uni [90] . In this case, components are defined by their phase space distribution f{p,x) 



11 



where p and x are respectively momentum and spacetime coordinates. Interactions are included as 
collision terms in the right hand side of the Boltzmann equation |9H [92| I93j . from which one can 
obtain energy-momentum and number conservation equations directly: 

^p^^ 



P^d^fiip, x) - Vp't^ ^ L[M = a{p, x) (34) 



nf^ = j dp a{p, x), dp^ -A_d'p6{E^ -f- m}) (35) 
Tru= ! dpp^C,{p,x) (36) 



where g is the number of internal degrees of freedom (e.g. spin) of species i. Conservation equations 
()35p and (|36p are obtained by using the following property of the Boltzmann operator L defined in 
see e.g. 



dp p^p^'^p^'^ . ..p'^"f{p, x) 



J dpp'''pf'K..p''"L[f{p,x)] 



(37) 



Collisional terms can be written by using cross-sections of interactions which can be determined 
separately from the quantum formulation of the model \91\ [39] . In the context of interacting dark 
energy models, the simplest examples of collisional terms are elastic scattering between dark matter 
and dark energy and slowly decay of dark matter with a small branching ratio to dark energy |36|^ . 
Note that we assume no interaction between dark energy and visible matter and radiation. For these 
interactions the collisional terms are: 



Cm{p,x) = -TmmmfmiP,x) - fm{P,x) J dp^ UiP<p,x) Ak{p,p^) (Jm<p{P,Pf) + 

dpm dp^ fm{Pm,x) f^{p^,x) Ak{prn,P^) '^""""^ '^'^'^^ (38) 
C^{p.,x) = Tm-mm j dpm fm{Pm,x) ^-^^P^'P^ _ f^{p^x) J dpm fm{Pm,x) Ak{p,Pm) 

dpm. dp^ fm{Pm,x) f^{p^,x) Ak{pm,P^) ""^ (39) 

Akipi,P2) = [{pi.p2f - m\ml]^ (40) 

where Tm is the total decay width of dark matter, M{pm,p) is the multiplicity of with momentum p 
in the decay remnants of dark matter particles with momentum pm-, and crm<piPm,P<p) is the total cross- 
section of interaction between dark matter and dark energy with momentum pm and p^, respectiveljlf]. 

The disadvantage of this approach is that it needs phase space distribution of components which is 
not always available, specially for dark energy. Moreover, the absence of an explicit description for the 
Lagrangian means that the total energy-momentum tensor needed for determining Einstein equations 
and metric evolution, can be obtained only by solving equation (|36p for all components. These equa- 
tions are differentio-integral and usually don't have analytical solution. Thus, in practice interacting 
models can be studied only numerically, otherwise one needs to consider some approximation. For 
instance, dark energy interaction with matter must be very weak. Thus, {Tj^^l <^ \ ^jZf^l- There- 
fore, we can neglect its contribution in the total energy-momentum tensor and Einstein eg nation J^. 



*In models where energy is transferred from dark energy to dark matter, the interaction must be nonhnear and very 
sophisticated such that a very hght quintessence field be able to produce massive dark matter particles. At present no 
fundamental description for such models is available. 

^Note that although dark energy is a condensate i.e. its particles have the same energy, presumably zero momentum, 
a general condensate state can contain very large number of particles in different energy levels, see [71] for more details. 

^''in some dark energy models such as early dark energy it is assumed that the density of dark energy at high redshifts 
is much larger, and only at low redshifts it is reduced. Although at redshifts relevant for dark energy surveys cosmology 
must be very close to ACDM, one must be aware that in many models of this type, the approximation of weak interaction 
can be applied only at low redshifts. It is also expected that these models leave a detectable signature on the CMB 
spectrum [57l[58] . 



12 



As for the integration of collision term in equations ()35p and ()36p . under some physically motivated 
assumptions they can be simplified and integrated. For instance, when dark matter is assumed to 
be an scalar, the expression for the scattering cross-section is very simple, see e.g. j36j . It is simply 
proportional to the coupling constant and delta functions for energy-momentum conservation. It is 
expected that the mass of quintessence field be very small, specially much smaller than the mass of 
dark matter. The momentums of both components are also expected to be small. In this case their 
distribution at large momentums is strongly suppressed, and the cross-section around the peak of 
distribution can be considered to be approximately constant. Under these simplifications, it is easy 
to see that scattering term in the right hand side of (j36|) is proportional to integrals of the form: 



/f TX^ XL f TV^ UjO Th' 

dpi dp2 Pi fl(,Pl,x) f2{p2,x) = "-1 / dp2 f2{P2,x) ^ ^ / dp2 P2 f2{P2,X^ _ 1 P 
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(41) 

nf^^, nr«^^ = ^ (42) 
\ni\ rrii rrii 

where nf and are number density and velocity of species i, respectively. Approximate expression 
for in ()42p is valid when the distribution in momentum space is concentrated around a peak. Using 
similar approximations the decay terms in the right hand side of (j36p can be also described as a 
function of velocity and number vectors. Finally, after grouping all the constant or approximately 
constant factors together, energy-momentum conservation equations for dark matter and dark energy 
can be written as: 



-^m ■,u ~ L^n^ -\- Amgn^^Uippn^ — Qm, ('^'^) 
-^(p ~ ^•fi'^m ~l~ ^ips^^Umpn^ = (44) 

where constants Li and Ais are decay width and scattering amplitude for species i. In the rest of 
this work we use these equations as an approximation for energy-momentum conservation equations 
irrespective of dark matter type (spin) and details of interaction between two dark components. They 
affect constants Li and Ais which are used as parameters. One can also add dark matter self- 
annihilation term to (03]). But, it is easy to show that self-annihilation is proportional to Ira^p. Thus, 
it is significant only in dense regions i.e. at small spatial scales such as the central region of dark 
matter halos, which are in nonlinear regime and are not studied in the present work. Here we only 
consider homogeneous and linear perturbations. Therefore, the effect of annihilation is negligible. We 
remind that equation (j43p is not restricted to cold dark matter and can be also used for relativistic 
matter, for instance neutrinos in early universe, or a hot component at low redshifts. 

Although in the rest of this work we consider the interaction terms described in this section, for what 
concerns the study of differences between modified gravity and interacting quintessence models, the 
formulation of anisotropics and discrimination methods explained in the next two sections can be 
applied to other choices of interactions. It is enough to find an interaction current similar to what we 
have found for decay and scattering above and add them to the right hand side of equations (|43p and 

dm). 



4 Cosmology and evolution of anisotropies 



In this section we first determine J-'i coefficients defined in Sec. [2] for both modified gravity and 
quintessence models according to interaction currents and energy-momentum conservation equations 
obtained in the previous section. Then, we consider the effect of interactions on the evolution of 
anisotropies, and describe how interaction parameters can be extracted from data. 
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4.1 Interaction coefficients in the Friedman equation 



4.1.1 Modified gravity 

Using the energy momentum conservation equation (j28p and the interaction current for modified 
gravity models, the scalar field equation and the evolution equation of the homogeneous matter density 
can be determined as the followings j67j : 

a' 
a 



^" + 2n^' + aV^(^) = a'^C{^) Y,{pi -^Pi)^ % = - (45) 



p[ + m{p^ + P^) = C{^W{pi - 3P,) i = m, 6, /i (46) 

where barred quantities are homogeneous components, in subscript means derivative with respect 
to ip. Note that here we have generalized the original calculation in [B7] by considering a (/9-dependent 
C{(p) coefficient in the right hand side of these equations to cover a larger class of modified gravity 
models, see e.g. [35] • For f(i?) models C = yMvrG/S [67]. Equations ()15|) and (jl6|) are coupled and 
an analytical solution can not be found without considering an explicitly V{^)- Therefore, to solve 
the equation for p, which is in fact the only directly observable quantity, we simply consider the right 
hand side of the equation as a time-dependent source. The solution of equation ()36]l can be written 
as: 

p,(z)=p,(zo)(l + ^)'('+"')e(i-3-')^('^), F{ip)^ [C{^)d^, i = m,b,h (47) 



where Wi = Pi/pi for all species except dark energy are assumed to be constant and are given in 
equation ()14p . Comparing this solution with ([3]) we find: 

T,{z) = e(i-3-0^{^W) « 1 + (1 - 3wi)F{^iz)) (48) 



In i(R) models C{ip) = y^AnG/S = C [67] is a constant, thus: 

Fi^) = C^iz). (49) 

Using transformation from Jourdan frame to Einstein frame ipiz) = ln(f/j(z) + 1)/2C [67], one can 
relate J-i(z) to Ir: 

T^{z) ^ 1 + in(f^ + 1) « (1 + f«) ^ (50) 

The approximate expression in (j50p is the same as equation (|14p. Note that in (j47p all constant 

coefficients including (1 + iR{zQ)) 2 are included in pi{zo). Apriori one can test the presence 
of a f(ii) modified gravity by measuring simultaneously Trn{z), Fh{z), and equation of state of dark 
energy from equation (jlSp . In fact in this equation if we neglect the last term that depends on the 
time derivative, the effective dark energy density becomes: 

3 fn , dlnfjR) 
^''^'^ SvrG 1 + ffi^ 6dR ^ ^ 

To be consistent with observations f(i?) cannot be a fast varying function of R. Therefore, the 
dominant term in ()5ip is the term proportional to R which makes the relation between p^e, R, and 
f{R) very simple. Other J^j's and evolution of corresponding densities have also known expressions, 
notably J-h{z) = 1. Therefore, apriori simultaneous fitting of these quantities can test f(i?) modified 
gravity models. More generally, in modified gravity models dark energy term in the Friedman equation 
is an effective contribution generated from non-conventional interaction between matter and gravity. 
Therefore, it is more correlated to matter than in ACDM or (interacting)-quintessence models. In 
the former apriori there is no correlation in the dark sector, and in the latter case the interaction can 
be very small and is only necessary for reducing fine-tunings and making the model more natural. 
Similar correlation tests can be performed for other modified gravity models too. Evidently, giving 
the small deviation of dark energy from a cosmological constant, the measurements and calculation 
of correlations are not trivial tasks. Furthermore, the discrimination must be cross-checked by using 
anisotropics for distinguishing between dark energy models, explained in Sec. 14.21 
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4.1.2 Interacting quintessence 



In the same way, we can determine Fi coefficients for (interacting)-quintessence using equation (|i3]) . 
We replace n'^ with approximation (j42|) and include 1/m factors in the L and As coefficients. After 
these simplifications, the evolution equation for the density of interacting quintessence models becomes: 

p'i + 37i{pi + Pi) = -Liapi + Asiapip^ (52) 

where i indicates any cold matter or relativistic species that interact with quintessence fielcJ^. A clear 
difference between interaction term in ()52p and (06]) is that the former does not explicitly depend on 
the scalar field, and therefore we do not need to know and solve a field equation similar to (|45p 0. 
The solution of this equation and corresponding J"j's are: 

p,{z) = p,(zo)(l + ^f^^+"'')exp(L,(r(z)-T(zO))+^, J ^^^"^(z) ) ^^^^ 
^iiz) = exp(^-Li(T(z) - t{zo)) + A,, j dzj^^^^^^ « 1 + Li(r(zo) - r(z)) + 

A,, rdz (54) 

where t{z) is the age of the Universe at redshift z. Note that even in absence of expansion, the density 
of dark matter at high redshifts is higher if Lj > 0. 

Along with consistency relation explained above for modified gravity models, explicit dependence 
of (|54p on measurable quantities p^{z) and H{z) apriori allows to discriminate between interacting 
quintessence and modified gravity models. Note that the prior knowledge about the evolution of these 
quantities are mandatory for distinguishing the underlying model and without such information one 
cannot single out any of these models. 



4.2 Matter perturbations in interacting dark energy cosmologies 

Although dark energy influences the evolution of perturbations mainly through quantities related to the 
homogeneous component - background cosmology - the study of anisotropies can be a powerful mean 
both for measuring the equation of state and for discriminating between candidate models. Standard 
candles, such as supernovae type la, allow direct measurements of distances, and thereby cosmological 
parameters. However, they are rare events, can deviate from being standard due to absorption or late 
detection [62], sub-types, and dependence of their light curve on other properties such as metallicity, 
mass, and magnetic field of progenitors [63] • Determination of dark energy properties from evolution 
of perturbations provides additional information and a mean for cross-check of the two methods. 

Matter perturbations in presence of an interacting dark energy [6l] and in f(i?) modified gravity 
models [19l [20l [671 EH] have been calculated by various authors, thus here we do not repeat them 
and simply use their results. Our main objective is to find and discuss features that can be used for 
discrimination between dark energy models. 

Considering only scalar perturbations, we define the first-order metric in conformal gauge as the 
following: 

ds^ = a^{r])[{l + 27/'(x))dr/2 - (1 - 2(j){:>c))6ijdx'dx^] (55) 

^^If species i has interaction with another component, for instance is scattered by another species, we can add a second 
scattering term to (|53|l . The best example is the scattering of photons or neutrinos by baryons. Here for the sake of 
simplicity we neglect such interactions which are not the main concern of this work. However, in a full formulation of 
the problem they should be considered, specially if they can mimic an interaction with dark energy. 

^^For f{R) modified gravity in which C is constant ip' in (|46|l can be replaced by an expression depending on density 
and pressure, and there is no need for solving field equation of the scalar field either. 



15 



As we mentioned in the Introduction, for modified gravity models we write evolution equations in 
Einstein frame. Thus, here only their interaction terms distinguish them from quintessence models. 

We use fluid description for both matter and dark energy. After linearizing energy-momentum con- 
servation equations and taking their Fourier transform with respect to spatial coordinates, evolution 
equations for density and velocity perturbations of matter component i and dark energy can be written 
as: 

6p[,) + 3n6p^i){l + + (1 + «;(i))p(,)(3</>' - ikjv^^^^) = ^Q(,)o (56) 

6p'^ + m5p^{l + C%) - (1 + w;^)p^(30' - ikjvip) = 6Q^o (58) 
((1 + W(i))p^v^j)' + 4H(1 + w^)p^v^j - ikjCl^5p^ - ikiU'-^j - ikj{l + w^)p^'ip = 6Q^j (59) 

where C^^^^ = dPii-^ / Spii-^ is the speed of sound for species i El, V(j_) is its velocity, and n'.^^. is its 
anisotropic shear. The perturbation of interaction current for modified gravity and quintessence 
models derived from ([30]) . ([I3]l and are as the foUowings: 



Modified gravity: 

^Q{i)o = P{i) 



Interacting quintessence: 



1 - 3w^,))C^{^)^'6ip + C{ip) ( (1 - 3«;(,)) V + (1 - SC^ ^))<^'%) 
—SQipo 

ikf{(p)pi^i){l - ?>wi^))5ip 



SQipj, i = m, b, h 



Pv^Pii) + P{i){^Pip + Pvi^) 



SQ{i)o = -aL(i)((5p(i) +/3(i)V') + a^s(i) 

^Qii)j = o.vi^i~^j{-Li^i~^p{^i) + As(^op(^i)p^) 

6Q^o = aL^{6p(^i) + ^jpli)) + aA^^ 5p^p(^) + p^5p(^i) - p(^i)6p^ ( - 

5Q^j = aV(^i)j{L^p(^^'^ + i = All matter interacting with (p 



1 + ^'^+^ 



(60) 

(61) 



(62) 
(63) 

(64) 
(65) 



In ([60|) . is the derivative of C{ip) with respect to if. To obtain these equations we have used the 
following definition and properties: 



Pn 



Sp^ + SP^ 
Pip + P,p 



d,^d-^ ~ {p^ + P^)-2 



{Pv + h 



(66) 
(67) 
(68) 

(69) 



Evidently, these equations are valid for both modified gravity and interacting quintessence. They 
are also highly coupled, thus it is impossible or very difficult to find an analytical solution for them. 
To complete evolution equations for modified gravity, we also need the evolution of 5ip. This can be 



To prevent confusion between spacetime indices and indices indicating the species, when there is a risk of confusion 
we put the latter inside brackets 
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obtained by expanding the field = ^ bip and using the covariant field equation, see e.g. [7T] : 

-^d^i^^gg^'^d.p) + V^{ip) = C{ip)Tm (70) 

V' + 2^ V + + a?^'^) + ^{^" + + 2— (^) - (fc2 _ ^ + y )5(^ = e{p)5T^ + c^f <5v? 

a a 

(71) 

As we mentioned in previous sections, solving these equations is not the main aim of present work. 
Our goal is to single out differences of these models that can be used for discriminating them from 
other models. For instance, in modified gravity models the perturbation of interaction current does 
not depend on the metric perturbations and (f). By contrast, in interacting dark energy the current 
perturbation depends on the metric perturbation and it is easy to that: 

termo^ ^ ^ .^2) 
term oc 5 

Because according to observations 6<p, <p' , and dp' are very small, in both models the terms proportional 
to 5i = 5 Pi/ Pi are dominant. In this case, it is easy to see that for modified gravity SQ(^i)Q oc C{p) 
and for interacting quintessence (5Q(j)o {—L^ + A^(j)p^). Although apriori these quantities evolve 
differently, both of them are expected to vary very slowly. Thus, it is not possible to distinguish them, 
specially in a model independent way. Other properties such as (f72]) cannot be used directly either. 
Nonetheless, they influence the growth rate oc SySi, density power spectrum, and density- velocity 
correlations, etc.. In the next section we discuss how these measurable quantities can be related 
to interaction current, and thereby allow to discriminate between modified gravity and quintessence 
models. 



Perturbation equations ()62p to ()65p depend on metric perturbations ip and (p, and their time derivatives. 
These quantities can be determined from Einstein equations for perturbations (see e.g. [95]): 



k^(l) + m{(l)' + nij) = ATTGa^'^5pi (73) 

i 

i 

cP" + nii;' + 2<p') + ( — - ^)V + ^ = -47rGa2 V 5P, (75) 

i 

k\ct> - ^) = -127rGa2 ^^(fc.A:' - ^4)% (^6) 



Note that in these equations the interaction energy is neglected. The reason is that we need T^^^^, 
which in the phenomenological description of interactions is not known. Nonetheless, its omission in 
equations (j73p to (j76p should not induce large errors because present observations show that any non- 
gravitational interaction between various constituents of the Universe - if any - must be very small, and 
therefore this approximation is justified. Metric perturbations ip and p cannot be directly observed, 
except through lensing. Otherwise, they can be extracted from these equations when density-density 
and density-velocity correlations, and induced anisotropic shear 11^^^^ are determined from LSS data. 

Although phenomenological interaction currents (j43p and (|44p are inspired from well understood scat- 
tering of particles, one cannot rule out other types of interaction. Even for these cases apriori one 
should be able to write equations similar to (|62p -(|65 p . and (|72p . The fact that the latter relations are 
independent of the strength of the coupling between dark energy and matter proves that finding a 
different proportionality between tp and 5 terms would be a clear signature of an unusual quintessence 
model, e.g. one with a non-minimal interaction with gravity. Evidently, such measurements are not 
easy. Nonetheless, with the huge amount of data expected from near future surveys and their bet- 
ter precision, more accurate measurements of parameters should be possible, and precision analysis 
necessary for detailed examination of dark energy models should be achievable. 
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5 Estimation of forecast precision for surveys 



In this section we first describe how in practice the background cosmology parameters defined in Sec. 
[2] are calculated. Their uncertainties determine how well a survey can discriminate between modified 
gravity and (interacting)-quintessence models, independent of the data type or observation method. 
Then, we calculate and parametrize the evolution equation of the growth rate of matter anisotropics 
and discuss its measurement uncertainty. As an example we make an order of magnitude estimate 
for the expected uncertainty of these quantities for the Euclid mission |43] . As we mentioned in the 
Introduction, a proper forecast needs detailed study of observational effects and uncertainties which 
is out of the scope of present work. 



5.1 Discriminating between a cosmological constant and other models 

As we discussed in Sec. 12.21 discrimination ability of surveys between a cosmological constant and 
a redshift dependent dark energy can be evaluated by using the function Q{z) defined in (|24|) . To 
calculate the quantity and its uncertainty, we need to know uncertainties of the estimation of effective 
background cosmological parameters. The function @ depends on ^^^j, w^^j{z), ^ej|) and w^^^^j{z), 
effective dark energy fractional density and equation of state dark energy determined, by fitting H{z) 
and A{z), respectively. By measuring H{z), from either supernovae or BAO data, one can determine 
w^J^j{z) and ^^^j relatively easily. On the other hand, measurements of w^J^^j{z) and ^^ejj are less 
straightforward, because one has to determine dH/dz, or equivalently dDA/dz and d'^DA/dz^ (see 
Appendix [SI for relation between these quantities). For this reason, the uncertainty of is dominated 
by uncertainties oiw^^^j{z) and ^^g//- FiriallY) coefficients FiS that present the evolution of equation of 
state of various constituents, are determined by fitting the deviation of H{z) from the null hypothesis 
of a ACDM cosmology. However, as we argued in Sec. [21 there are strong degeneracies between J^j's 
and 7(z) which can be resolved only with using other types of data, in particular matter anisotropies, 
see Sec. 15.21 for more details. 



As an example, we estimate the uncertainty of G for the Euclid mission. For the parametrization 
Weff{z) = 'Wp + ■Waz/{1 + z), according to the Euclid-Red Book [9l|, the standard deviation for these 
coefficients are expected to be a^p ~ 0.015 and ~ 0.15 for Euclid data alone, and a^p ~ 0.007 
and a^a ~ 0.035 for Euclid+Planck data. No forecast for the expected uncertainty of dH/dz is yet 
available. For this reason, we simply use error propagation rules to determine a rough estimation 
for (Jiiijiiiz from available forecasts. We approximate dH/dz with its definition as a difference ratio: 
dH/dz ~ AH/Az, then we use the general uncertainty propagation rule to a function of n variables 

/(xi, . . . , Xfi)'- 

^ ^^^^ 

i,j=l,...,n ■' 

where Cij is the covariancG matrix for random variables xi, • • • , ^n- 

Assuming ou/H ~ 1%, neg- 
ligible error for and Ti ~ 1, the dominant source of error in dH/dz from wi^z). Because the 
coefficients of derivatives with respect to these parameters in (j77p is roughly of the order of one, 
we estimate a i / {dH / dz) ~ 10 — 15%. Functions A{z) and B{z) are related to dH/dz, see 
(jlOSp . and when the uncertainties of H, density fractions Jlj's and redshift z are much smaller, 
(7^{A) /w^^'^l{z) ~ '^qW /^i// ~ '^o"s ~ ^dH/dz ~ 10% around optimal redshift of z ~ 0.5. Measurement 
precisions of J^j's also are of the order of precision of i.e. ajjTi ~ a j / {dH / dz) ~ 10 — 15%. 

Evidently, uncertainties obtained here are very rough estimations. The aim of these exercises is just 
to show what level of error we expect from near future surveys. A proper prediction needs detailed 
simulation of measurements and data analysing methods, instrumental effects, and systematic and 
statistical errors. They need a dedicated study that we leave to a future work. 

Finally we want to make a remark about the redshift dependence of w{z), which in the literature 
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is usually parametrized [5^ . In Appendix |A] we show that for the same value of w at two different 
redshifts, different parametrizations lead to very different evolution for A{z). Inversely, if we estimate 
w{z) from the measurement of A(z), parametrization of w{z) can lead to very different evolution for 
this function, despite employment of the same data for A{z). Therefore, we must estimate w at each 
redshift without parametrizing it. As for the estimation of uncertainties, for instance from the Fisher 
matrix, they can be determined from the set of {w{z),'y{z), z} at every redshift bin rather than from 
a parametrization, see Appendix [B] 



5.2 Discrimination between modified gravity and interacting quintessence models 

If we observe a non-zero 0, then we must use the power spectrum and growth rate of perturbations to 
investigate the nature and origin of deviation from a cosmological constant. The comparison between 
evolution equation of modified gravity and interacting quintessence models in Sec. 14.21 showed that 
their interaction currents are very different, and thereby the evolution of matter anisotropies and 
dark energy density in these models are not the same. In fact, if we could decompose the interaction 
current to terms proportional to scalar metric perturbations and matter density fluctuations, it were 
possible to distinguish between these models. However, in practice measured quantities are matter 
power spectrum and its growth rate i{z,k) defined as: 

f(,,t)^'-^^:^, (78) 

a ma rtdm Om{z = 0,k) 

The function f{z,k) is usually extracted from the power spectrum using a model [961 E3 EH]) for 
instance a power-law for the primordial spectrum including its modification by Kaiser effect \99\ [TOOl 
lion 1102] and redshift distortion due to the velocity dispersion |103j . 

To obtain the evolution equation of f(z, k), we replace potentials V ^ind (p by expressions depending 
only on 5m = Spm/pm and 9^. = 'i-f^jV'j^m)' Assuming a negligible anisotropic shear at z < 0(1) 
which concerns galaxy surveys, scalar metric perturbations - gravitational potentials - ip and (j) can be 
determined from Einstein equations (|73| ) - (f75| )P^ 

<P = ^P = ^^^(Sm + 3{l + w„,)'^)+A^P (79) 



AttG / 



Ai; = -^i^dp^-m6ip{p^ + P^)2 j (80) 
AirGprnnf, ,^ k'^ .nOn 

= — 1^ ■ ■ 

= -nAtP + AirGaH^ip^ + P^)^ (82) 



6m + {S + ^)(1 + ^m)^] + A0' (81) 



Note that in ()79p and ()8ip we have separated terms which vanish for ACDM model and written them 
as A-i/^ and A(p'. As observations show that dark energy behaves very similar to a cosmological constant 
- at least for z < 0{1), both these quantities are expected to be very small. It is why we write them 
as a variation of -0 and 0'. For future use it is also better to redefine them as followings: 



1 



_ 6p^ _n{p^ + p^)26(f 

eo = ^— , ei = = (ocij 

Pm Pm 

M = ^^{e,-^e,) (84) 
= p (eo-(3 + ^)eil (85) 



^^In this section for the sake of simplicity of notation we consider that Fi's factors for species are included in uii's, i.e. 
(1 + zf'<^Ti is redefined as (1 + z)3Ti(^) and Wi is obtained from Q using this redefined 7i. Therefore, for interacting 
dark energy models Wm is nonzero and in general depends on redshift. 
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After replacing (j)' and -0 in ([56|) and ([57|) with ([8T]) and ([79]) respectively, evolution equation of matter 
and velocity perturbations can be written as: 



5'^ + ^ + m{{l + CDSm + (1 + Wm] 
Pm I 



<5m+ (3+ + m 



+ 



(eo - (3 + ^) [ + (1 + wm)e„, = 6Qmo (86) 

= ihSQ\^) (87) 

where 6Qmo and ikiSQ^^^^ are interaction currents and 0^ = SvrGa^/Jm/S^^- Moreover, in present and 
near future wide area surveys such as DES and Euclid the value oi'H/ck <^ 10. For instance, for Euclid 
T-L/ck < 0.01. Therefore, we can neglect terms proportional to H/k. Under these approximations, 
evolution equations of density and velocity become: 

Modified gravity: 

5'^ + 3n{C]^ - Wm)Sm + (1 + Wm)Om = 



1 



8.G «^^i+C(^)^ ^ ) 2n,il + w,f ^'' 

0'm + - ^f^6m - ^Ti^S^ + 60 - 3e,) = _ ^k'(^-^^rn)nm c^^^,^ (89) 

1 + 2 {8ttG{1 + w^)n^)2 

Interacting quintessence: 

aAsrn^Q (90) 

9'm + - ^^Sm " ^n^{Sm + 60 - 3ei) = -^^(-L,„ + AsmP^)aem (91) 

l + Wm 2 l+Wm 

Now that we have the evolution equations for 6m and 9m, we can determine the evolution of growth 
rate. The procedure for calculating di(z, k)/dz is straightforward. We replace 9m in (|89p and ()9ip with 
its value obtained from (j88p and (|9ip . respectively for modified gravity and interacting quintessence 
models. Then, we replace 6'm with its value from equation ()78p . The final equation has the following 
general form: 

ru + f (H' + n^) + f^H^ + 3(^2^ - wm){n' + m^) + 3(^72^ _ ^ 

^ nm{l + u;™)^^^ + A;2C2^ + E^m + + E2U + E^U^ + Ei = {) (92) 

Coefficients Eq, Ei, E2, E3, £'4 depend on z and k, and have the following values for the two models 
discussed here: 



^^Note that the speed of light c = 1 is assumed in metric (|55|l . and therefore it does not expHcitly appear in our 
calculations. 
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Modified gravity: 



El 
E2 



C{Lp)a{p^ + P^) 2 y3{Cgj^ - Wm) + 1 - ^Wr, 



3(1 + Wm){Cg^ - Wm)C{ip)a{p^ + Pf)^ - C^{ip){l - 3Wr^ 



E3 



30m(l + Wmf ( ep 3ei" 



+ 3(C7; 



(93) 
(94) 

(95) 
(96) 



I 



C^{ip){\ - 3w 

m I 



n 



+ C{ip)a{p^ + P^y2 



i{l-3wm){l + C^Jpmeo\' 



2(1 + w^)p^ 



+ 



(3C{(p)a{p, 



^ + P^)2{Ci^ - W.m) - C((^)a(l - 3w.m){p^ +P^)^ X 



'l-?,Wm){l+ClJ-p 



2(1 + w^)p^6„ 



Wr, 



Interacting quintessence: 



Eq 

El 
E2 

E3 
Ea 





SWm.aiCg^ - Wm){-Lm. + AgmPip) + ^smaPm(l + 3w„ 

3J^m(l + u!mf ( eo 3ei 

2 \ 5m 

- f e'o 2aeo. 4 - 



(97) 

(98) 
(99) 
(100) 

(101) 

(102) 



In the calculation of ()92p - ()102p we have neglected terms proportional to Ti/ck. 



For ACDM model Ei = 0, i = 0, • • • , 4. For a non-interacting quintessence model all Ei coefficients 
are zero except E3. A notable difference between modified gravity and interacting quintessence models 
is the coefficient Ei which is strictly zero for interacting dark energy models and nonzero for modified 
gravity that leaves an additional scale dependent signature on the evolution of matter anisotropics. 
The other explicitly scale dependent term is common for all models and is expected to be very small 
because it is proportional to the square of sound speed which is very small for cold matter. In addition, 
in contrast to the rest of Ei coefficients, Ei and E3 are dimensionless. Evidently, the contribution of 
Eik^ term with respect to other terms in equation (j92|) increases for larger k, i.e. at short distances. 
But, the effect of nonlinearities, i.e. mode coupling also increases at large k, see e.g. |104j . They 
can imitate interactions and lead to misinterpretation of data. For this reason, it is suggested that 
observation of galaxy clusters is a good discriminator between dark energy models |105[ llOd] . because 
clusters are still close to linear regime, but have relatively large k. 

Discriminating power of a survey can be estimated by the precision of Ei and £'3 measurements. 
However, one expects some degeneracies when equation (j92p is fitted to determine Ei^s. Moreover, 
in galaxy surveys, f and f (or more exactly di jdz) are determined from the measurement of power 
spectrum from galaxy distribution, and'H and from the BAO effect on the spectrum. Thus, these 
measurements are not completely independent. An independent measurement of 7^ and 'H! e.g. using 
supernovae will help to reduce degeneracies and error propagation from measured quantities to the 
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estimation of Ei^s. The relation between T-L' and B{z) defined in (jllOp shows the logical connection of 
parametrization of homogeneous component - background cosmology - and evolution of fluctuations, 
specially in what concerns discrimination between dark energy models. In fact, anisotropics depend on 
the equation of state of matter, which in the context of interacting dark energy models, is modified by 
its interaction with dark energy. Thus, their independent measurements optimize their employment 
in distinguishing between various models. 

Although apriori df /dz can be determined directly from data by differentiating f , usually due to shot 
noise the errors would be very large unless we extensively rebin the data. However, rebinning smears 
the redshift-dependence, which is the most important information for discriminating between models. 
Another approach is to solve equation (j92p analytically. It does not have an analytical solution for the 
general case, but as we show in Appendix[Dl when Wm and are approximated by constant values, 
and cosmology is matter, radiation or cosmological constant dominated, i.e. up to desired precision 
only one component determines its evolution, an approximate solution can be found. At present epoch 
where matter and dark energy have comparable contributions, coefficients in ()92p even for ACDM 
vary with redshift. Nonetheless, their variation arrives very quickly to saturation. Therefore, the true 
solution is not very different from the approximate analytical one under the explained conditions, and 
it is possible to determine perturbations around the analytical solution by linearizing equation (|92p . 
see Appendix [D] for more details. 

A rough estimation of the uncertainties of E'j's measured by Euclid can be performed in the same 
manner as what is presented in Sec. 15.11 for Q and J^i's. It is expected that growth rate f can be 
reconstructed from Euclid+Planck data with an uncertainty cif/f < 3% [M]. Considering equation 
([92|) and estimation of uncertainty of Ti' obtained in Sec. 15. H the uncertainty of (Jf /f must be 
~ 10% — 15%. This limits our ability to distinguish between a ACDM model where Ei = Wm = 0, 
and quintessence or interacting dark energy models where these quantities are not zero. Considering 
the linear equation obtained in Appendix |D] from expansion of f around its solution for ACDM, the 
total uncertainty of deviation from this model is roughly the same as what is obtained for f, i.e. 
~ 10% — 15%. But, the uncertainty in the estimation of each Ei is expected to be larger because 
of the degeneracy of these parameters. Evidently, determination of f and f at multiple redshifts 
should help somehow reduce degeneracies and improve discrimination between models. More precise 
estimations as well as the estimation of the effect of nonlinearities and the optimal choice of scale 
range need detail simulation of surveys. We leave these tasks for future works. 

5.3 Interpretation and comparison with other parametrizations 

It would be useful to have a better insight on the physical meaning of the parameters defined in the 
previous section, and to compare them with what is used in the literature for parametrizing dark 
energy models. 

We begin with eg and ei defined in ()83p . Their definitions show that the former depends only on dark 
energy density anisotropy and the latter only on the peculiar velocity of dark energy field, i.e. on its 
kinematics, see (|124p . They follow each other closely and approach zero when the field approaches its 
minimum value. However, their exponent close to the minimum depends on the interaction. There- 
fore, their measurements give us information about the potential and interactions of the scalar field. 
Moreover, the difference in the dependence of evolution equation of anisotropics and growth factor to 
eo and ei shows that only by separation of kinematics and dynamics of dark energy - scalar field - it 
would be possible to distinguish between modified gravity and other scalar field models. 

The deviation of gravity potentials (p and ip from their value in ACDM A?/^ is the quantity which 
can be measured directly from lensing data |107j . For this reason various authors have used A^/J to 
parametrize the deviation from ACDM [TBI [771 [751 [7S] . However, equations ([7^ and ([SU]) show that 
although Atp ^ is by definition a signature of deviation from ACDM, in contrast to claims in the 
literature, it is not necessarily the signature of a modified gravity model because quintessence models. 
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both interacting and non-interacting, also induce l^ip 7^ 0. This is also another manifestation of 
the difference between kinematics and dynamical effects of interacting dark energy models described 
above. 

Because we have used Einstein frame for both quintessence and modified gravity models, in absence of 
an anisotropic shear (p = ip even in non-ACDM models. In linear approximation gravitational lensing 
effect depends on the total potential ^ = (j) + ip (see e.g. |lU7j for a review). Therefore, in Einstein 
frame 

^ = 2<P = 24^ = $ACDM + ^ACDM = {^^ + 3(1 + Wm)^^ = (103) 

In the notation of [76] $ = 2E$acdm, thus: 



^ACDM k^^m 

The other quantity which affects the evolution of lensing and directly depends on cosmology is the 
growth factor of matter anisotropies which determines the evolution of defined in (|103p . This 
quantity can be obtained from integration of growth rate f defined in ()78p and is usually parametrized 
as flJn- For ACDM 7 w 0.55 [108] . In this respect there is no difference between our formulation 
and what is used in the literature. Evidently, this simple parametrization cannot distinguish between 
various dark energy models. By contrast, the more sophisticated decomposition proposed in Sec. 15.21 
is able to distinguish between quintessence and modified gravity. Note that in Jordan frame there 
are two other parameters: rj = {ijj — (p) / cj) and Q = (p/(j)ACDM- The parameter S = Q(l + rj/2), thus 
it is not independent. In Einstein frame rj = 1 unless there is an anisotropic shear. At first sight it 
seems that there is less information in Einstein frame about modified gravity than in Jordan frame. 
However, one should notice that in Einstein frame the fundamental parameters are eo and ei and 
other quantities such as A^l' and f can be explained as a function of these parameters. Therefore, the 
amount of information in Einstein and Jordan frame about modified gravity - if it is what we call 
dark energy - is the same. The advantage of formulation in Einstein frame and definition of eo and ei 
is that they can be used for both major categories of models. Moreover, they have explicit physical 
interpretations that can be easily related to the underlying model of dark energy. 

More recently based on an original work by C. Skordis |109j . two groups [llOUlll] have suggested new 
parametrizations which are basically only for discriminating modified gravity models from ACDM. 
Both groups use the following approximate description for the Einstein equation: 

G^u = R^u - \9f.uR = SttGT^, + U^, (105) 

The tensor U^u is called Energy-momentum tensor of dark energy |109| . and originally its definition 
has been for formulation of all modifications of the Einstein theory of gravity. In [llOj this tensor is ex- 
panded with respect to potentials ip and (j), and coefficients of this expansion are used for parametrizing 
the underlying modified gravity model. 

Note that equation ()105p is at all scales an approximation because the right hand side is explicitly 
proportional to the Newton coupling constant. Considering f{R) models which are the simplest 
modification of the Einstein theory of gravity, in contrast to (|105p , the coupling to matter is modified 
in both frames, see equations ([8|)- ()10p for Jourdan frame, and the formulation of f{R) model in 
Einstein frame in |67] . In fact, in Einstein frame the modification is explicit in the energy- momentum 
conservation equation. This means that if a deviation from ACDM is observed, it would be very 
difficult to verify the consistency of the model at short distances because the deviation of coupling 
from Newton constant G is put by hand to zero. Moreover, this formulation and parametrization by 
definition does not help to detect interaction between dark energy and matter, because it depends 
only on the total variation of metric potentials. In addition, in this formulation Ufj^i, is assumed to 
be a conserved component, which as we discussed in Sec. [3l is not consistent because in contrast e.g. 



23 



to perturbative quantum field theories, we never measure the free component. Furthermore, equation 
(|105p has exactly the same form for quintessence models, thus in this framework it is not possible to 
discriminate between this class and modified gravity models without knowing the underlying model 
in detail. 

The formulation in uses a Lagrangian formalism with quadratic and higher order deviations from 
the Einstein theory of gravity. The energy-momentum tensor of dark energy Ufj,i, is obtained be using 
variational methods from this Lagrangian. It is a function of g^i, or the set {g^u, ^, dfj_<p} when 
the dark-(energy) sector includes also a scalar field. Then, they use 3+1 spacetime decomposition, 
thus all coefficients of the above expansion depend only on time, and apply variational methods to 
determine perturbations SU^i, around an arbitrary background. Their formulation is technically and 
theoretically interesting, specially for studying various modified gravity models, but there is neither a 
model independent parametrization for dark energy nor for observables. 

6 Outline 

We have parametrized the interaction between dark energy and matter for modified gravity and 
interacting quintessence models as modifications of the evolution of matter and radiation background 
and perturbations densities, and the equation of state of dark energy. We have showed that when the 
interaction is ignored in the data analysis, the effective value of parameters are not the same if we 
calculate them from Friedman equation or from a function proportional to the derivative with respect 
to redshift of total mean energy density of the Universe. We have also defined a single quantity that 
evaluates the strength of the interaction. Its observational uncertainty can be used to estimate the 
discriminating power of a cosmological survey. 

We have obtained a phenomenological description for the interaction current in the context of inter- 
acting quintessence models motivated by particle physics. Based on these results, we have suggested 
to distinguish between modified gravity and (interacting)-quintessence dark energy models of non- 
gravitational origin by the way they modify energy-momentum conservation equation. If the inter- 
action current is proportional to the trace of the energy-momentum tensor of matter, we classify the 
model as modified gravity, otherwise, as (interacting)-quintessence and its variants, such as K-essence, 
quintom, cosmon, etc. 

We have determined the modification of evolution equation of density and velocity perturbations in 
the context of modified gravity and interacting quintessence models discussed above, and used them 
to obtain a parametrized description of evolution equation of the growth factor that can be used for 
both these models as well as a simple ACDM model, which has been considered as the null hypothesis 
in our discussions. The difference between the value of these parameters can distinguish between 
aforementioned models. We have also obtained order of magnitude estimations for uncertainties on 
these quantities measured with the Euclid mission. A better forecast for these uncertainties needs 
simulations of the survey and the data analysis that we have left to future works. 
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A Properties of A{z) 



One of the principle aims of LSS surveys is the measurement of Hubble constant H{z), angular diameter 
distance Da, and luminosity distance Dl, mainly by measuring Baryon Acoustic Oscillations (BAO) 
which play the role of a reference distance scale [IB]. The maximum effect of BAO on the power 
spectrum is at redshift ~ 0.3 [l6j. However, as we mentioned in the Introduction a direct determination 
of 7(2:) from Hubble constant. Da, or Dl when z — )• is not possible. In fact, using equation ([1]) and 
the definition of angular diameter distance, it is easy to see that: 



In 
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Inside + ^l{z)\og(l + z) (106) 



At small redshifts the last term on the r.h.s. of ()106p which contains 7(2) approaches zero, and the 
effect of the latter becomes negligibly small. Now, consider the following quantities: 
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where H{z) = a/a is the expansion rate of the Universe and p{z) is the total density at redshift z. 
It is clear that A{z) is proportional to the deviation of dark energy from a cosmological constant 
at any redshift including z = 0. In addition, its sign determines whether dark energy has normal 
or phantom-like equation of state at a given redshift. It can be shown [37j that when dw/dz <^ 
3w{z){w{z) + 1)/(1 + z), the sign of dA/dz is opposite to the sign of w{z) + 1. This condition is 
satisfied at low redshifts - see examples of models in Fig. [TJ It means that A{z) is a concave or 
convex function of redshift, respectively for positive or negative w{z) + 1. Observations show that the 
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contribution of ilk and 0,h at low redshifts is much smaller than the uncertainty of f^m- The function 
dA/dz does not depend on ^rn- Thus, the uncertainty on the value of Q^m can shift the value of A{z) 
but it does not change its slope and its shape i.e. its concavity or convexity will be preserved. 

The function B{z) can be easily related to directly measurable quantities: 

, _J_(dDi _ J2j^\ _ d^Di 
Biz) = = ^ = ^+"^ '^^ i±£^ ^ fllO) 

dz 

°' = <' + -''^°X flW 

(112) 

or equivalently with respect to normalized angular distance: 

Biz) - -(^^ + a + ^)4g^) ni3) 

Ho r dz Di 

''^ = TT-zJolm = OTW ^'''^ 

Note that these equations are written for a flat universe, but can be easily extended to the cases where 



B About Fisher matrix for equation of state of dark energy 

Fisher matrix evaluates the sensitivity - information content - of a measured quantity to variables 
and parameters that define the underlying model [IT]. Under special conditions, e.g. Gaussianity of 
distributions. Fisher matrix can be related to the covariance matrix of measurements. In LSS surveys 
the main measured quantity is the power spectrum of matter density anisotropies. Application of 
Fisher matrix to CMB jJH] and galaxy surveys [961 [971 [HE] is well studied and widely used. In 
what concerns the measurement of dark energy density, its variation, and its equation of state from 
galaxy surveys, one has to extract H{z) and Da{z) either from BAO [112[ I113[ I114j or by fitting the 
complete power spectrum [51]. Fisher matrix for 2-dimensional power spectrum is determined by Seo 
& Eisenstein [961 [971 [98] with H{z) and Da{z) as parameters. A transformation from these quantities 
to coefficients of a parametrized equation of state, for instance w{z) = wq + Waz/{1 + z) allow to 
determine the covariant matrix for the measurement of wq and Wa [51| . 

Although apriori the value of these quantities can be determined at any redshift, in practice the limited 
volume and deepness of surveys allow to determine the power spectrum at the average redshift of the 
survey or for some bins of redshift in the case of large deep surveys. In the latter case, the estimation 
of w{z) as a function of redshift depends strongly on its parametrization. Fig. [T]shows the plot of A{z) 
for examples in which w is measured at two redshifts. It is evident that this quantity and thereby the 
underlying dark energy models depend strongly on the parametrization of notably when systematic 
and statistical errors are added. 

Simpson and Peacock [115| use {wQ,Wa^^K,^k^^mh'^ ,^bh'^ ^n^^ Ag, (5,^' ,ap] as independent param- 
eters for estimating cosmological parameters from the measurement of the galaxy power spectrum. 
Here Wa = —dw/dz, /3{z) = fiz)/hiz) where i{z) is the growth rate of scalar fluctuations and b{z) is 
the linear bias, and 7' is the parameter that define an approximate parametrization for f(z) « f^m(-z) 
for ACDM [108] . It can be also shown that in what concerns the determination of the Fisher matrix 
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Figure 1: A{z) as a function of redshift. To see how well A{z) can distinguish between various models 
and how systematic and statistical errors as well as parametrization affect the reconstructed model, we 
consider 3 parametrizations as written on the plot above. Note that parametrizations for the plot in the 
center and on the right are equivalent up to a redefinition of coefficients wq and wi . We first consider a 
given value for w{z) at z = and z = 3, determine corresponding coefficients WQi and wu where index 
i is for initial. Then to simulate systematic errors we plot the following models: wq = —1 + \wQi + 1|, 
wi = wii (dotted line), wq = —1 — \woi + 1|, wi = —wu (dot-dash), wq = —1 + {woi + 1|, wi = —wu 
(dashed) and wq = wqi and wi = wu (full line). Colored vertical bars present statistical errors. The 
uncertainty of az is lo"^(2=o) = 0-01 (top row) and lcr^(2=o) = 0.05 (bottom row) at z = and evolves 
with redshift as (Ja{z) = (ta[z = 0)(1 + z)^. It seems to be possible to distinguish between normal and 
phantom dark energy models easily, if uncertainties are limited to few percents. Evidently, achieving 
such a precision is challenging even for space missions such as Euclid. 
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for dark energy, w{z) and dw/dz alone lead to a singularity Qj. 

In place of parametrizing 'w{z), we suggest to use w{z), 7(2) and z to determine the Fisher matrix 
for dark energy parameters. It can be easily shown that Fisher matrix becomes singular if the first 
two quantities are considered [51], because w{z) and 7(2) are not independent - if one knows w{z), 
then 7(2;) can be determined from ([2]). This problem does not arise when w is parametrized because 
expansion parameters are explicitly independent. The relationship of w{z) and 7(2;) is very similar 
to the relation between H{z) and Da{z). Fisher matrix for {H{z), Da{z), z} set of parameters is 
calculated in |51j . Using this formulation, a parameter transformation gives the Fisher matrix for 
{w{z),'y{z), z}. Relation between Fisher matrices with 2 sets of parameters pi and qm is 



- \ - dqm dqn , ^ 



For the parameter-sets discussed above, the components of the Jacobian matrix are: 
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Alternatively, one of w^z) or 7(2) parameters can be replaced by A{z') = Qdei^iz) + 1)(1 + z)^^"'~^\ 
In fact, it is preferable to replace j{z) with A{z), because at low redshifts the 7(z)-dependent term 
has very small effect on the evolution H{z) and Da- By contrast, the deviation of A{z) from its value 
in ACDM model is maximized for z — t- 0, see Fig. [TJ 



C Fluid description of a scalar field 



Energy momentum tensor of a scalar field is: 



= --g''''g'"'dpipd,ip + g^^'Vi^) + d^^d-^ (122) 



Using definition ([55]) of a perfect fluid, the density and pressure are defined as: 

= u^u.T!^^ = ^d'^^d^if + V{v), ^ \d^^d^,^ - y{^) (123) 

is the velocity vector and u^u^ = 1. It is easy to verify that with above definitions for and P^p: 



1 



(124) 



^^For the sake of simplicity in the discussion of Fisher matrix here, we neglect other cosmological parameters, i.e 
we assume that dark energy parameters can be factorized from other quantities. In practice, one has to consider a 
single matrix Fisher matrix containing all parameters. Thus there would be one single covariant matrix that includes 
correlation of all uncertainties. 
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D Solution of evolution equation of growth rate 



For ACDM cosmology, Ei = 0, 



0, . . . , 4. We also consider Wr, 



dividing equation (j92p by % , the evolution equation of growth rate becomes: 



^ + f(^ + i) + f^ + ^o. 







After changing the variable from r] to In a, this equation changes to: 



0, 



In 



By integrating the Friedman equation for flat ACDM one obtains: 



n 



8 = ^ 



drj 



3x 



p/2 



-3x 



0. In this case after 
(125) 

(126) 

(127) 
(128) 



For z = 0, i5 = — 1 and for z ^ oo, £ = ^h{a) — ^rn{o)- To be able to solve ()126p analytically we 
must assume is a constant. This is a good approximation if we are interested only on a small range 
of redshifts. Under this assumption, the solution of ()126p can be obtained by integration: 



^kCDM{z) 



-{£ + 1 



^) + {£ + l + 



ai = V(f + 1)2 - 60, 



(129) 



1 - (1 + z)°i 

For 1 < £ < ^/QQrn — 1) IS imaginary and according to this approximation solution f(z) 

has an oscillating component. A simple attempt to make (|129p more precise is to take into account 
that £ depends on redshift. 

To obtain an approximate solution for interacting dark energy models parametrized by coefficients 
-Ej = 0, i = 0, ... ,4 in (|92p . under the assumption that these corrections are small, we can linearize 
this equation around (acdm- Note that in general it is expected that in interacting dark energy 
models Wm and are not zero. Therefore, we add also their contribution to the linearized model: 



f = ^ACDM + Af 

'n' 



(130) 



Af' + 



n 

3T~L I C?, 



■y2 

■'sm 



+ 2 1 



ACDM 



Wr. 



+ —Wm{2 + Wr, 



Af + 3(C: 



2 

sm 



Ea 



+ {Ctm + Ei)—+E2 + E3n + -^ = (131) 



Solution of this linearized equation is straightforward and can be formally written as the following: 



Af(z) 



n 



(l + z)(l + 3(C7L -t/;„)) 

(l + z)(l + 3(C2 



exp 



dz 



1 + 

m 



dz 



l + z 

Wm)) 



-L 9 f 

+ ^ lACDM 



2 

sm 



wm)^ + m[{c^ 



1 + Z)H2 

(( 



Wm) + 



exp 



dz ( Ei 



n 



l + z 

Wm{2 + Wr, 



+ 2 ^ACDM 



fe2 

n 



(Csm + El) + E2 + 



Ezn + 



n 



(132) 



Determination of integrals in (|132p needs details of redshift dependence of coefficients EiS which is 
model dependent. Nonetheless, they depend on the scalar field which must vary very slowly with 
redshift. Therefore, at zero order, they can be considered as constant. Although even with this 
simplification it is difficult to determine (|132p analytically, a numerical determination allow to write 
it as an expansion with respect to Ei coefficient. This expansion would be suitable for compression 
with data and determination of E^. 



34 



